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SUMMARY 

This paper considers the estimation of longitudinal aerodynamic coefficients from data recorded dur- 
ing flight tests of a powered- lift STOL aircraft. First, a comparison is made between the coefficient 
values determined by the regression and quasi linearization identification techniques from records taken 
during elevator pulse maneuvers. The results show that for these tests the regression method provides 
less scatter in coefficient estimates and provides better correlation with the predicted values. Special 
techniques are then developed which allow identification of the coefficients from records taken during 
landing maneuvers in which the aircraft encounters turbulence while flying in ground effect. Flight test 
results are presented to illustrate the effects of air turbulence and ground proximity on the estimated 
coefficient values . 

NOMENCLATURE 

pitching acceleration, rad/sec 2 T 

a x acceleration measured along X-axis, g units u 

a 2 acceleration measured along Z-axis, g units V 

c mean aerodynamic chord, m w 

C aerodynamic coefficient x 

g acceleration of gravity, m/sec 2 a 

h height -above- ground- level, m <5 

I yy inertia about the Y-axis e 

K constant parameter p 

M aircraft weight a 

q pitching rate, rad/sec » 

Q dynamic pressure 

S aircraft wing area, m 2 


1. INTRODUCTION 

NASA is conducting a rather broad research program on powered- lift concepts for future use with jet 
STOL transport aircraft. As part of this program a C-8A Buffalo aircraft has been modified with an aug- 
mented jet-flap system (ref. 1). This aircraft has been undergoing flight tests to determine the in-flight 
aerodynamic performance and handling qualities. In support of this program a study has been made to 
evaluate the use of parameter identification techriques in determining the aerodynamic coefficient values 
from the recorded flight data. 

Several identification methods are available from previous studies (refs. 2-10) to identify the air- 
craft parameters from the records taken where the aircraft is excited only by elevator Inputs in calm air. 
These previous methods, however, are generally unable to treat the problems associated with identification 
of the aircraft parameters during landing maneuvers where there are significant external disturbances due 
to the air turbulence and ground proximity. 

In this investigation two different parameter identification techniques have been applied to data 
recorded during pulse- type maneuvers where the aircraft dynamics are excited by elevator inputs. This 
paper will review the accuracy in determining the coefficient values using these different identification 
techniques. Special techniques are then applied to data recorded during landing maneuvers where the air- 
craft it excited by the combination of air turbulence, ground proximity, and the pilot's normal control 
actions. This paper reviews the development of these special techniques and presents results which Illus- 
trate the effects of air turbulence and ground proximity on the estimated coefficient values. 

The intent of the paper is to present the genesis of each of the problems and the identification 
algorithms used in the problem solution along with a discussion of some of the more important findings. 


thrust term 

velocity along X-axis, m/sec 

total velocity, m/sec 

velocity along Z-axis, m/sec 

vector of state variables 

angle-of -attack, rad 

elevator deflection, rad 

pitch angle, rad 

atmospheric density 

standard deviation (rms) 

free air value, out-of-ground effect 

estimated value 
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2, AIRCRAFT AND INSTRUMENTATION SYSTEM 

The results in this paper ware obtained fro* flight tot data recorded during tot maneuvers with an 
augmented jet-flap STOL research aircraft (ref. 1). Thi* vahicla (fig. 1) is a high-wing STOl aircraft 
powarad by two turbofan anginas mounted in nacallas located undar the wing. The relatively cold flow fro* 
the front fans is ducted to augmentor jet flaps. The angina exhaust is directed through nozzles, one on 
each side of the nacallas, to provide vectored propulsive lift. 

The flight test instrumentation included a nose boon with a pitot-static systea and vanes, body- 
mounted accaleroaet ers and rate gyros, vertical gyros, position transducers on the control surfaces, pres- 
sure and temperature transducers to measure the propulsive characteristics, and a radar altimeter to 
measure height- above- ground- level. The vane-measured angle-of-attack, a, has been corrected to account 
for angular rates and for upwash (as a function of hei ght- above- ground) . The pitching angular accelera- 
tion, a., has been derived from the pitch rate signal. The linear accelerations, a z and a x , have been 
obtained from the body mounted accelerometer signals and corrected (to the aircraft center-of -gravity) to 
account for angular accelerations. The flight data were obtained with an airborne digital recorder and 
then processed at discrete points, 10 points/ sec, on a ground based digital computer. 

3. COMPARISON OF ESTIMATION TECHNIQUES 

This section will review some estimation resuits for standard pulse-type maneuvers in which the air- 
craft is relatively free fro* turbulence effects and is above ground proximity effects. Emphasis will be 
to compare results fro* the different identification techniques and to gain some understanding of their 
relative accuracy in estimating the values for the aerodynamic coefficients. 

Several previous studies (refs. 2-10) have compared different identification algorithms for estimat- 
ing aircraft parameters and have found that the results may depend on the technique used. These identifi- 
cation techniques generally fall into two categories: equation error and output error. With noise in the 

measured aircraft states, the equation error technique can produce biased estimates of the coefficient 
values (refs. 2-4). The output error technique can reduce the bias error; however, it is affected by 
modeling errors and also may produce the larger standard deviations in the estimated coefficient values 
(ref. 10). This paper will compare results of both a regress on technique (equation error) and a quasi- 
linearization technique (output error). 

3.1 Identification algorithms 

The non-linear equations used to mathematically model the aircraft longitudinal forces and pitching 
moment were taken as: 


a x * (QS/M) [e Xo ♦ C Xfl a ♦ C X{ t ♦ C^(q£/2V)] ♦ T x (1) 
i x - (QS/M) [Cj q ♦ C la a ♦ C Z( t ♦ C lq (qc/2V) ♦ C^cJ ♦ T, (2) 
i, • (QSc/lyy)^ ♦ C^a ♦ C„ 4 « ♦ q^(q£/2V)] ♦ T, (s) 


Coefficient terms are included which account for variations in the aircraft angle-of-attack, a; elevator 
deflection, 6; and pitch rate, q. This model also includes a C z term due to the powered-lift function, 
Cj, (Cj ■ thrust of cold air/QS). Using this model the unknown coefficient values have been determined by 
the regression (also called equations of motion, or least squares) and the quasi linearization (also called 
modified Newton-Raphson) parameter identification methods. (Reference 10 outlines the details of these 
techniques as used for the results in this report.) 

Regression is a relatively simple technique which determines the coefficient values that minimize the 
least squares difference between the time histories for each of the measured accelerations, a x , a z , and 
a*, and the corresponding model outputs, 4 X , & t , and a*. The coefficient values are determined in three 
independent solutions, eqs. (l)-(3), using the well-known matrix inversion procedure (ref. 7). 

Quasi linearization, in contrast to the regression method, integrates the following kinematic equa- 
tions to obtain estimated time histories of the aircraft states. 


u ■ ,(», - (in 5) - qw ♦ K u , 

u(o) « u 0 

(4) 

W • g(»j ♦ co$ 8) ♦ qu ♦ £* , 

w(o) - w 0 

(S) 

♦ 

7 

• 

« ■ O’ 

q(o) » q. 

(6) 

J • q ♦ Kg * 

»(o) ■ »o 

(7) 


This technique determines the coefficient values (and bias terms) that minimize the weighted least squares 
difference between the time histories of the measured variables, a x , a z , a*» u, w, q, and 9, and their 
corresponding estimated values. With this technique, Initial estimates for the unknown parameter values 
are made (e.g., from the regression results) and then the estimates are successively improved in an 
iterative manner, using the quasilinearization algorithm (refs. 2, 3 and 6). 
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One primary difference between these two nethods is that with the regression method the variables, q, 
a, V, and Q, in eqs. (l)-(3) are taken as the measured values, whereas, with quasi linear! ta t ion thes e 
variables in eqs. (1)«(3) are represented by the estimated values 4 , q, a ■ tan* 1 (w/u), $ ■ /u* ♦ 5*, and 
$ • pV z /2. 

3.2 Comparison of estimated and measured time histories 

Figures 2 and 3 present a coaparison of measured time histories with those computed using the two 
identification methods. Figure 2 presents the regression results and fig. 3 presents the quasilineariza- 
tion results. Values for the rms difference between the measured and estimated data are listed in table 1. 
Figures 2 and 3 illustrate that the estimated time histories generally fall within the scatter of the 
measured data. As shown in table 1, the rms fit to the pitching acceleration, a*, is about the same for 
both methods; however, the regression method provides as much as a 30% better fit to the measured linear 
accelerations, a, and a 2 . 

3.3 Comparison of coefficient values 

The coefficient values determined by the two techniques are presented in fig. 4. Also shown (dotted 
lines) are the corresponding values which have been predicted from other independent sources, such as 
steady-state flight tests, wind tunnel tests, and theory (refs. 1, 11-13). 

In general, the more important coefficients such as C z , C x , C* , and , are in agreement both 

o a 6 q 

between the two methods and with the predicted values. The standard deviations (e.g., run-to-run scatter) 
of these estimated parameters are also relatively small. 

Other coefficients, such as C 2 , C 2 , C x , and C x , show somewhat more scatter. The inability of 

6 q 6 q 

either technique to estimate these terms accurately is probably because the influence of these terms on 
the aircraft forces is small. And also, there is a strong dependency between the elevator deflection, 6, 
and the pitch rate, q. Previous studies (e.g., refs. 7, 14-16) have also noted the large standard devia- 
tion associated with estimating these terms. 

For almost all of the coefficients, the regression values have less run- to-run scatter and agree 
better with the predicted values. A majority of the regression values (with the exception noted above) 
are within about ±10% of the predicted values. 

3.4 Discussion of identification techniques 

The results presented show that the regression method provides better results than the quasi lineariza- 
tion method. For instance, the regression method provides a better fit to the measured accelerations, 
less scatter in the estimated coefficient values, and better agreement with the predicted values. 

Any errors to be expected with the regression method depend, to a large extent, on the amount of 
measurement noise. Any noise in the measurement of the variables, q, a, V, or Q, could cause bias errors 
with the regression method. Although the amount of noise cannot be determined with certainty, the 
recorded data (e.g., fig. 3b) show very little of what may be termed white or near white measurement noise 
(e.g., there is a low noise- to- signal ratio). Apparently, for the flight test situations considered in 
this study, there are no large amounts of measurement noise that could cause significant errors with the 
regression method. 

The errors to be expected with the quasi linearization method are related to inaccuracies in the esti- 
mates of the variables, q, a, V, and Q (fig. 3b). In particular, any modeling enrors (e.g., neglect of 
higher- order aerodynamic terms and cross -coup ling from the lateral -directional mode) will cause inaccu- 
racies in these estimated states. Also, the quasi linearization technique usually has larger standard 
deviations in the estimated coefficient values because all of the coefficients are determined within one 
dependent set of equations, eqs. (l)-(7); whereas using the regression method, the coefficients are 
determined with three independent equations, eqs. (l)-(3). 

For this particular application, regression appears to be the better method to use in obtaining the 
coefficient values. This should not imply that in other situations (i.e., where there may be larger 
amounts of measurement noise, or where all the states are not directly measured) regression would be the 
better method to use. Experience has shown that it is good practice to consider both methods utilizing, 
wherever possible, the advantages of each method. 

4. PARAMETER IDENTIFICATION IN TURBULENCE 

One ef the problems in parameter identification during landing maneuvers is to account for the air 
turbulence which is usually present at low altitudes. Most of the previous studies have considered air- 
craft parameter identification in the absence of turbulence (refs. 2-10), or have made simplifying 
assumptions about the noise spectrum of turbulence and its interaction with the airframe (refs. 17-18). 

In this paper, a state estimation technique (ref. 19) is used to measure the time history of the turbu- 
lence gust disturbances. This measured turbulence is then treated as a forcing function in the aero- 
dynamic equations. This technique makes no assumptions about the turbulence noise characteristics and 
further allows an examination of the manner in which the turbulence interacts with the airframe. 

4.1 State estimation 

The inertial velocities and position of the aircraft have been estimated by a solution of the follow- 
ing kinematic equations; 



T 


ft • q(a x - sin §) - 

qw ♦ 

V 

u(o) ■ K u 

0 

(8) 

w ■ g(a z ♦ cos 0) ♦ 

qu ♦ 

H * 

w(0) . K* 

0 

(9) 

0 ■ q 

♦ 

V 

0(0) - fc e 

0 

(10) 

A « u sin 0 - w cos 

6 ♦ 

\ • 

h(o) - i ho 

(11) 


where the unknown constant (K) terms are determined by quasi linearization. This application of quasi* 
linearization requires no mathematical model of the aerodynamics; rather, the direct measurements of the 
cce 1 era t ions, a z and a x , and the pitching rate, q, are used in a manner similar to that in a strapped - 
I down inertial system (see ref. 19 for further details and the formulation including lateral motions). 

For the landing approach maneuvers in this study, state estimation provides smoothing of the measured 
states, (h, d), along with the estimates of the inertial velocities, (u, w), and the inertial angle-of- 
attack, ■ tan* 1 (w/u). Figure 5 presents some of the estimated states along with the measured data for 
a representative landing approach maneuver. The upper portion of the figure illustrates good correlation 
between the radio altimeter measurement and the estimated height -above-the-runway. The lower portion of 
the figure compares the estimated inertial angle of-attack, d^, and the airflow (vane-measured) angle-of- 
utack, a-. For these representative landing app*-oach maneuvers there seems to be a large random fluc- 
tuation or the airflow vane. The difference between the airflow and inertial angle-of-attack provides a 
measure of turbulence acting on the airframe. 

Note : An examination has been made to determine possible errors in the airflow angle-of-attack 

measured ”"by the vane. As noted previously, the airflow measurement, a a , includes corrections for angular 
rate and upwash (as a function of height -above- ground) . It has been found that for flight maneuvers out 
of turbulence there is excellent agreement between the airflow measurement, a a , and the estimate, d^. 

4.2 Interaction of turbulence with the airframe 

A necessary consideration in parameter estimation for STOL aircraft, traveling at low speeds, has 
been to determine the manner in which this measured turbulence interacts with the airframe. The turbu- 
lence as measured by the vane located forward of the aircraft will not immediately interact with the major 
aerodynamic surfaces. A first approximation for this delayed interaction is to use a time shift, At, to 
account for the time it takes for the measured gusts to travel from the vane until they strike the major 
aerodynamic surfaces. 

Noting that the total angle-of-attack at any time consists of both the gust and inertial components, 
we have 

i Measured turbulence Inertial 

shifted by At angle-of-attack 

a(t) * a gust^ * 4 «i(%) ( l2 ) 

I 

where the turbulence gust component is obtained as the difference between the measured airflow angle-of- 
i attack and the inertial angle-of-attack at the time, t - At. 

a gust(t - At) « a a (t - At) - - At) (13) 

\ Figure 6 illustrates the effect of the tine shift. At, on the rms fit errors, o a ^, o t , and o a • foT 

l a typical segment of a landing approach maneuver. As shown, there appears to be a different value of time 

; shift. At, which will provide a minimum rms fit error to each of the measured terms, a x , a z , and a^. 

These values of time shift appear to be reasonable from aerodynamic considerations. The fit error for the 
linear forces, o. and o. , are minimized if the measured turbulence is delayed by the amount of time 
l x 

required for the gusts to travel from the vane to near the aircraft aerodynamic center (At * 0.4 sec at 
V ■ 36 m/sec). The fit error for the moment term, o^, however, is minimized using the time required for 

the turbulence to reach the stabilizer (At * 0.7 sec at V » 36 m/sec). 

, The relative amounts of rms fit error reduction, with the time delay, also appear reasonable. The 

linear z force is strongly affected by angle-of-attack gusts and, as shown by using the appropriate time 
shift, the rms error, o a # is reduced by about 30%. The moment term and the linear x force are influ- 
enced less with a reduction of about 10% in and 5% in by the appropriate choice of tine shifts. 

j A further indication of the importance of time shift becomes apparent in fig. 7 where the effect of 

At on the estimate values for the coefficients, C, , fi* , and C. , is shown. Without a time shift (at 

1 q * 

; ! At * 0) the estimated values are much different than predicted. However, using an appropriate time shift 

* I (At * 0.4 sec) these terms are near their predicted value. 

; t 

4.3 Discussion of turbulence effects 

The appropriate value of time shift is related to the ratio, length/speed. For the linear forces, a z 
and a x , the time shift can be taken approximately as: 


at « distance from vy to ilrcraft A.C. 
forward airspeed 
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For the pitching moment the time shift is approximately: 

_ distance fTom vane to stabilizer 
At Forwirl affspeeJ 

With large STOL aircraft flying at low speeds the appropriate time delay will be in the order of seconds. 
For small aircraft at high speeds , however, the time delay may be quite small. 

Previous parameter identification studies, which have included turbulence effects, apparently did not 
find a requirement to time -correlate the vane-measured turbulence. These previous studies (refs. 17-18) 
have considered smaller aircraft at higher speeds where the inclusion of the time shift may not be so 
critical. However, as shown by the results in this paper, the time -dependent interaction of turbulence on 
the airframe can affect significantly the estimated coefficient values and, therefore, should be considered 
in each application. 

One additional note is that turbulence may, in fact, aid in the identification of some of the param- 
eters. This is because turbulence acts as another forcing function in addition to the usual control input 
forcing function. The results from this study indicate that some of the aerodynamic coefficients may be 
determined more accurately from maneuvers in turbulent air (e.g., the aircraft is excited by both gusts 
and elevator inputs) as compared with maneuvers inclear air (excited by only elevator inputs) . For 

instance, as noted previously, the terms C z and C 2 are highly correlated and difficult to determine 

q 6 

accurately using elevator pulse maneuvers (fig. 4). However, the estimated values in turbulence are 
generally found to be near their predicted values (fig. 7). 

5. PARAMETER IDENTIFICATION OF GROUND EFFECTS 

Ground proximity effects are of concern with STOL aircraft because wind tunnel tests and theory have 
predicted significant changes (both static and dynamic) in the aerodynamic flow field for such high- lift 
aircraft near the ground (refs. 20-23). These effects on the aerodynamic coefficient values have not yet 
been determined by accurate in-flight measurements from landing maneuvers. This section reviews a pre- 
liminary application of parameter identification to determine the changes in the aerodynamic coefficients 
due to ground proximity. Parameter identification has been used in two ways. First, it has been used to 
determine the gross changes in the aerodynamic coefficients due to ground effect. Second, it has been 
used in the development of a mathematical model which indicates the amount of change in the aerodynamic 
coefficients as a function of height -above -ground, angle- of -attack, etc. 

Representative maneuvers, which have been used to analyze the ground effects, are presented in fig. 8. 
In each of these runs the pilot controlled the aircraft near a constant angle-of- attack. Maneuvers are 
shown at different levels of angle-of-attack for different nozzle angle settings (i.e., different levels 
of aerodynamic and propulsive lift). 

5. 1 Gross effects of ground proximity 

The gToss effects of ground proximity on the aerodynamic coefficients can be isolated as follows: 


% - C L - [C,J 

(14) 

ACp ■ Co - fC„J 
G 00 

(15) 

^Mq * c m - 

(16) 


where the terms, AC^, ACp G * and ACm g » represent the gross changes due to ground effect; the terms, Cl, 


Cp, and C^, are the measured aerodynamic coefficients, 

C L ■ [-(»z “ T z ) cos a ♦ (a x - T x ) sin a] (M/QS) (17) 

Cp ■ [-(a* - T z ) sin a - (a x - T x ) cos a] (M/QS) (18) 

C M - (»» - T«) (lyy/QSc) (19) 

and the terms, [C^], (fip ]# ] are the predicted coefficient values (total sums) derived from 

parameter identification out of ground effect (as discussed previously). 

ICiJ * \ * \ “ * % 5 • • • • 

[Cp 1 ■ Cp ♦ Cp a ♦ C D 6 * . . . (21) 

* o a o 

(Cm ] ■ Cg| ♦ Cg a ♦ C§ 6 ♦ . # • (22) 

• o a o 


Figure 9 present t representative results showing the gross changes in aerodynamic coefficients as a 
function of he ight-abcve- ground level. An examination of the data presented in thi» figure provides 
insight into some of the variables which influence t*' changes in the aerodynamic coefficients and also 
indicates the type of terms which must be included ie mathematical model for ground effect. 
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First, the magnitude of the ground effects generally vary in an exponential manner as the aircraft 
nears the ground. This type of variation with height is similar to that noted in most previous studies of 
ground effect. 


Second, the ground effects vary from run to run depending upon the aircraft operating conditions. 

For run 1 (shown by circle symbols) there is a more positive change in lift and a more negative change in 
drag ms compared with run 2. These differences apparently account for the greater increase in flight path 
angle and speed near the ground in run 1 as compared with run 2 (fig. 8). For run 1 the ground effect 
appears "buoyant" enough to cause the aircraft to float up away from the ground; whereas with run 2, the 
ground effect appears less buoyant, and the aircraft continues to descend to the ground. 

Third, fig. 9 shows that the magnitude of the ground effects is somewhat different for descent and 
ascent (shown by the arrows). This apparent "dynamic" grour. effect i: illustrated in more detail with 

fig. 10 where the time history of ACt is presented for the portion of run 1 where the aircraft descends 

G 

and ascends above ground level. As shown, there is a rather abrupt loss in lift associated with the change 
from the descending to ascending flight path. This decrease occurs after the passage of the minimum alti- 
tude point. Apparently, the effect of the ground plane on the flow field is time -dependent. Near the 
ground the flow field is effectively straightened, causing a lift loss (see sketch in fig. 10). Such a 
lift loss, with a time lag, has been predicted from previous small scale dynamic tests (ref. 20); however, 
it had not yet been verified from actual flight test data* 


5.2 Mathematical model for ground effect 


An examination of the data presented in fig. 9 (and similar data f^rci other runs) gives insight into 
the form of equations required to model mathematically the changes in t»e aerodynamic coefficients due to 
ground proximity. A preliminary mathematical model which is being evaluated is of the form 



iC °G 


e -h/K h 


(23) 


Kl 

where the tern e' h ^ K *' represents the exponential variation of the ground effect with height; [£] repre- 
sents a matrix (3xn) of unknown constant coefficients; and x represents a time varying vector (nxl) of 
state variables which influence the amount of change in the aerodynamic coefficients due to ground effect 
(e.g., angle-of-attack, rate of descent, etc.). 


The parameter which has been found to have the most significant effect on the rms fit error is scale 
height parameter, K^. Figure 11 illustrates the relative rms values for Cl, C&, and Cm as a function of 
Kh* As shown in fig. 11, the best fit is obtained, for all three coefficients, with a scale height param- 
eter of * 4.5 meters (15 ft). 

Using the values obtained by parameter identification we can see how each of the variables (e.g., h, 

а, etc.) affect the aerodynamic coefficients. As an example, fig. 12 presents the estimated aerodynamic 
coefficients as a function of angle-of-attack both in and out of ground effect. Ground proximity is shown 
to cause (1) a slight increase in C, at low angles -of- attack along with a slight decrease in the lift 
curves slope. Cl , (2) a reduction of about 30% in C D , and (3) a significant shift in the moment, Cm, 
with an increate°in the static stability, -C^ . 

The trends, due to ground proximity, found in this flight test study are in general agreement with 
results found in a wind tunnel study using a similar powered-lift STOL configuration. That is* the wind 
tunnel tests also show similar changes in lift and lift curve slope, a decrease in drag, and similar shift 
in moment with increased static stability. However, the magnitude of the changes are somewhat different 
in the flight tests as compared with the wind tunnel tests. Figure 13 compares the changes due to ground 
effect, AC^, ac D g * AC M G * as obtained from flight and wind tunnel tests. In comparing these data the 

Haight above the ground level has been normalized with respect to the chord length; also ACl and ACq 

G G 

are normalized with respect to their free-air values. As shown, the changes in lift and moment are in 
general agreement with the wind tunnel, however, the decrease in drag determined in the flight test is 
about three times greater than the decrease in drag determined in the wind tunnel. Some differences were 
to be expected between the flight and wind tunnel results because, in the wind tunnel, the angle-of-flow 
between the ground plane and airframe is not the same as in actual flight; and also, in the wind tunnel 
there is a boundary layer on the ground plane (for fixed planes), again not the same as in actual flight. 
Because of the difficulties of accurately duplicating the ground proximity effects (both static and 
dynamic) from wind tunnel tests alone, it would appear that parsmeter identification, as used in this 
study, can be an important tool in the analysis of ground effects for future vehicles. 

б. CONCLUDING REMARKS 


This paper has reviewed some recent flight experience in the identification of longitudinal aero- 
dynamic coefficients for a powered- lift STOL aircraft. Comparisons were made between results obtained by 
the regression and quasi linearization identification techniquas. Also, special techniques were presented 
for the identi fleet ion of aerodynamic coefficient! when the aircraft encounters air turbulence and 
ground proximity. 

This study shows that for the data analyzed in this investigation tha regrasslon method provides 
better results than the quasi linearization method. The regression method provides a better fit to the 



measured accelerations, lest scatter in the estimated coefficient values, and betteT agreement with the 
predicted values. 
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The technique for estimating parameters in turbulence involves the use of state estimation, combined 
with airflow (i.e., vane) measurements, to determine the time history of the gust disturbances. The 
results show that the measured turbulence must be time-correlated to account for interaction of the gusts 
along the airframe. Using this technique, the results indicate that some of the aerodynamic coefficients 
may be determined more accurately from maneuvers in turbulent air (e.g., the aircraft is excited by both 
gusts and elevator inputs) as compared with maneuvers in clear air (excited by only elevator inputs) . 

In the estimation of ground proximity effects parameter identification has been used in two ways. 

First, it ,ias been used to determine the gross changes in the aerodynamic coefficients due to ground 
effect, and second, it has been used in the development of a mathematical model for ground effect. The 
results show that ground proximity causes a slight increase in lift, a moderate decrease in drag, and a 
significant change in pitching moment. 

This review illustrates that there are some differences between the results obtained by the various 
identification methods, but of more importance, is a determination of the form of the aerodynamic equa- 
tions (i.e., number and type of nonlinear and time-dependent terms) required to model mathematically the 
aircraft and its interaction with external forces. For this study of STOL aircraft, during landing 
maneuvers 5n turbulence, the primary consideration has been to define the form of the mathematical models. 
Future work appears warranted to investigate the problems of developing the most accurate mathematical 
models for advanced STOL and V/STOL aircraft. The development of these mathematical models requires an 
analysis of the recorded flight data along with an understanding of those physical processes which may 
affect the vehicle dynamics. 
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